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Abstract. For amorphous carbon several laboratory ex- 
tinction data are available, which show quite a wide range 
of differences due to the structural complexity of this ma- 
terial. We have calculated self-consistent dynamic models 
of circumstellar dust-shells around carbon-rich asymptotic 
giant branch stars, based on a number of these data sets. 
The structure and the wind properties of the dynamical 
models are directly influenced by the different types of 
amorphous carbon. In our test models the mass loss is 
not severely dependent on the difference in the optical 
properties of the dust, but the influence on the degree of 
condensation and the final outflow velocity is considerable. 
Furthermore, the spectral energy distributions and colours 
resulting from the different data show a much wider spread 
than the variations within the models due to the variabil- 
ity of the star. Silicon carbide was also considered in the 
radiative transfer calculations to test its influence on the 
spectral energy distribution. 

Key words: hydrodynamics - radiative transfer - Stars: 
atmospheres - Stars: carbon - Stars: AGB and post-AGB 



1. Introduction 

Asymptotic giant branch (AGB) stars show large ampli- 
tude pulsations with periods of about 100 to 1000 days. 
The pulsation creates strong shock waves in the stellar 
atmosphere, causing a levitation of the outer layers. This 
cool and relatively dense environment provides favourable 
conditions for the formation of molecules and dust grains. 
Dust grains play an important role for the heavy mass 
loss, which influences the further evolution of the star. 

Condensation and evaporation of dust in envelopes of 
pulsating stars must be treated as a time-dependent pro- 
cess since the time scales for condensation and evaporation 



are comparable to variations of the thermodynamic con- 
ditions in the stellar envelope. The radiation pressure on 
newly formed dust grains can enhance or even create shock 
waves leading to more or less pronounced discrete dust 
shells in the expanding circumstellar flow (e.g. Fleischer 
et al. |i"99l] p9l| Hofner et al. p^95[ Hofner & Dorfi |l997|). 
Since a significant part of the dust grains transferred to 
interstellar space are produced in the atmosphere of these 



old luminous stars (Sedlmayr 1994) an understanding of 
the nature of mass loss of these long-period variables is 
crucial for the general understanding of dust in space. 

Modelling circumstellar envelopes requires knowledge 
of the absorption properties of the different types of grains 
over the relevant part of the electromagnetic spectrum. 
For this the optical properties of the corresponding dust 
material are needed. Amorphous carbon is a very good 
candidate as the most common type of carbon grains 
present in circumstellar envelopes, since the far-infrared 
data of late-type stars show a spectral index as expected 
for a very disordered tw o-dim ensional material like amor- 
phous carbon (Huffman 1988 ). 

Silicon carbide (SiC) grains seem to be another im- 
portant component of the dust in circumstellar envelopes. 
While amorphous carbon could explain the continuum 
emission, SiC particles could be responsible for the 11.3 
/im band observed in many C-rich objects. 

In this paper, we present self-consistent dynamical 
models of circumstellar dust shells calculated with selected 
laboratory amorphous carbon data. Based on these mo- 
dels we have performed radiative transfer calculations for 
pure amorphous carbon and in some cases also including 
SiC dust. In Sect. [| the used amorphous carbon data are 
described. The influence on the model structure is de- 
scribed in Sect.|] and the resulting spectral appearance 
is discussed in Sect. ^. 
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2. Optical properties of dust 

The two sets of quantities that are used to describe opti- 
cal properties of solids are the real and imaginary parts of 
the complex refractive index m = n + ik and the real and 
imaginary parts of the complex dielectric function (or rela- 
tive permittivity) e = e 1 + ie". These two sets of quantities 
are not independent, the complex dielectric function e is 
related to the complex refractive index, m, by e' = n 2 — k 2 
and e" = 2nk, when the material is assumed to be non- 
magnetic (/i = hq). Reflection and transmission by bulk 
media are best described using the complex refractive in- 
dex, m, whereas absorption and scattering by particles 
which are small compared with the wavelength are best 
described by the complex dielectric function, e. 

The problem of evaluating the expected spectral de- 
pendence of extinction for a given grain model (i.e. as- 
sumed composition and size distribution) is essentially 
that of evaluating the extinction efficiency Q ex t- It is the 
sum of corresponding quantities for absorption and scat- 
tering; Q cxt = Qabs + Qsca- These efficiencies are func- 
tions of two quantities; a dimension- less size parameter 
x = 2Tta/\ (where a = the grain radius and A = the 
wavelength) and a composition parameter, the complex 
refractive index m of the grain material. Q a bs and Q sca 
can therefore be calculated from the complex refractive 
index using Mie theory for any assumed grain model. The 
resulting values of total extinction can be compared with 
observational data. 

A limit case within the Mie theory is the Rayleigh ap- 
proximation for spherical particles. This approximation is 
valid when the grains are small compared to the wave- 
length, x — 2ira/X << 1 and in the limit of zero phase 
shift in the particle {\m\x « 1). In the Rayleigh approx- 
imation the extinction by a sphere in vacuum is given as: 



comparisons of laboratory measurements with astronomi- 
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2.1. Measuring methods of optical properties 



A proper application of Mie theory to experimental data 
requires that the samples are prepared such that the parti- 
cles are quite small (usually sub-micrometer), well isolated 
from one another, and that the total mass of particles is 
accurately known. 

In order to obtain single isolated homogeneous par- 
ticles, the grains are often dispersed in a solid matrix. 
Small quantities of sample are mixed throughly with the 
powdered matrix material e.g. KBr or Csl. The matrix 
is pressed into a pellet which has a bulk transparency in 
the desired wavelength region. Some of the problems with 
this technique are that there is a tendency for the sample 
to clump along the outside rim of the large matrix grains 
and that the introduction of a matrix, which has a refrac- 
tive index different from vacuum, might influence the band 
shape and profile. This matrix effect can be a problem for 



cal spectra (Papoular et al. 199S; Mutschke et al. 1999). 

By measuring the sample on a substrate (e.g. quartz, 
KBr, Si or NaCl) using e.g. an infrared microscope, the 
matrix effect can nearly be avoided since the sample is 
almost fully surrounded by a gas (e.g. air, Ar or He). But 
the amount of material in the microscopic aperture re- 
mains unknown, which is an important disadvantage of 
this method. Therefore, these measurements are not quan- 
titative but they reveal the shape of the spectrum nearly 
without a matrix effect (Mutschke et al. 1999] ). 

A major problem of both methods is clustering of the 
grain samples either during the production of the partic- 
les or within the matrix or on the substrate. Clustering 
can cause a dramatic difference in the optical properties 



(Huffman 1988). A way to avoid this problem is to per- 
form the optical measurements on a polished bulk sample. 
For the determination of both n and k two or more mea- 
surements on bulk samples are required. This might be 
done either by a transmission and a reflection measure- 
ment, or by two reflectance measurements determinations 
at different angles or with different polarisations. Since the 
real part, n, of the refractive index, m, is determined by 
the phase velocity and the imaginary part, k, by the ab- 
sorption, a transmission measurement easily fixes k. The 
Kramers-Kronig relations can be applied in order to ob- 
tain the optical constants for grain measurements. The 
real part of the refractive index can be expressed as an in- 
tegral of the imaginary part (see e.g. Bohren & Huffman 



1983) 



3. Carbon grains 

While carbon is expected to constitute a major fraction 
of the circumstellar dust in carbon stars, its exact form 
is still unclear. Carbon has the unique property that the 
atoms can form three different types of bonds through sp 1 , 
sp 2 (graphite) and sp 3 (diamond) hybridization. 

A number of observations of late-type stars contradict 
the presence of graphite as the dominant dust type (e.g. 
Campbell et al. 1976| ; Sopka et al. 1985: Martin & Rogers 



|1987| ; Giirtler et al. ^996j). The far-infrared (FIR) data 



of late-type stars generally show a dust emissivity law of 
Q(X) ~ \~P with a spectral index of (5 « 1. While graphite 
grains have a FIR emission proportional to A~ 2 (Draine 
& Lee |1984 ), a A" 1 behaviour can be expected in a very 
disordered two-dimensional material like amorphous ma- 



terial (Huffman 1988; Jiiger et al. 1998) 



Amorphous carbon grains therefore seem to be a very 
good candidate as the common type of carbon grains 
present in circumstellar envelopes. Another possibility 
could be small diamond grains. Presolar diamond grains 
have been identified from primitive (unaltered) meteorites 
(carbonaceous chondrites) and are the most abundant 
(500 ppm) of the presolar grains discovered so far (Lewis 
et al. 1987] ). At least 3% of the total amount of carbon 
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present at the formation of the Solar System was in the 



form of diamonds (Huss & Lewis 1994). The place of origin 
of the presolar diamonds is still unknown, but since they 
can only have formed under reducing conditions J0rgensen 
( |1988D has suggested C-rich AGB stars as the place of for- 
mation of the majority of the presolar diamond grains. 



It has been suggested by Kriiger et al. (1996) that the 



surface growth processes on carbonaceous seed particles in 
circumstellar envelopes will take place at sp 3 bonded car- 
bon atoms rather than at sp 2 bonded ones, which suggests 
that the grain material formed in circumstellar envelopes 
could be amorphous-diamond like carbon. Presolar dia- 
monds extracted from meteorites have a median grain size 
of about 2 nm (Fraundorf et al. 1989 ), meaning that each 



diamond contains a few hundred to a few thousand carbon 
atoms. The presolar diamonds therefore seem to actually 
consist of a mixture of diamond (core) and hydrogenate 
amorphous carbon (surface) having about 0.46 the volume 
fraction of pure diamond (Bernatowicz et al. 1990). 

Several spectra of presolar diamonds from various me- 
teorites have been published (Lewis et al. 
et al. |1995|; Mutschke et al. |1995| Hill et al. 



198E 



Colangeli 
1997| Andersen 

et al. 19981; Braatz et al., submitted to Meteorit. Planet. 



Sci.) and even though a number of artifacts tends to be 
present in all the spectra, the general trend is that the 
presolar diamonds have an absorption coefficient that is 
twice that of pure diamond and a factor of a hundred less 
t han t he "diamond-like" amorphous carbon of Jager et al. 
(|1998|) . 

There exists a wide variety of possible amorphous car- 
bon grain types, which fall in between the categories 
"diamond-like" and "graphite-like" amorphous carbon. 
We have calculated dynamical models using various labo- 
ratory data of amorphous carbon to determine the possi- 
ble influence of these different grain types on the structure 
and the wind properties of C-rich AGB star models. 

3.1. Laboratory measurements of amorphous carbon 



Laboratory conditions are far from the actual space con- 
ditions where grains are produced or processed, but ex- 
periments in which physical and chemical parameters are 
controlled and monitored do give the option of selecting 
materials which may match the astronomical observations. 
When choosing which amorphous carbon data to use one 
is faced with the fact that due to the various processes 
used in the sample preparation, differences often appear 
between the measurements of various authors. Another 
major problem is that the optical properties of amorphous 
carbon are most often obtained by different techniques 
in different wavelength regions. Extinction measurements 
of sub-micron-sized particles is the most common tech- 
nique in the infrared. In the visible and ultraviolet, reflecti- 
vity and transmission measurements are often obtained on 
bulk samples. 
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Fig. 1. Calculated extinction efficiency of amorpho us car - 
bon from optical c onstan ts published by Maron ( 1990 ), 
Roule au fc Martin (|l99l|) , Pr eibisc h et al. (|l993|) , Zubko 
et al. ( 1996 ) and Jager et al. ( 1998 ), see Table |l| for anno- 
tations. The extinction efficiency factors were calculated 
from the optical constants (n and fc) in the Rayleigh ap- 
proximation for spheres. 



Bussoletti et al. ( 1987a ) have determined the extinc- 
tion efficiencies for various types of sub-micron amorphous 
carbon particles and spectroscopically analysed them in 
the wavelength range 1000 A - 300 /im. In their paper 
they present an updated version of the data already pub- 



lished from 2000 A to 40 (Borghesi et al. |1983| , |l985a| ) 
and new data obtained in the UV/vis (1000 - 3000 A) and 
in the FIR (35 - 300 /xm). The sub-micron amorphous 
carbon grains were obtained by means of two methods: 

(1) striking an arc between two amorphous carbon elec- 
trodes in a controlled Ar atmosphere at different pressures 
(samples AC1, AC2 and AC3; where the numbers refer to 
different accumulation distances from the arc discharge); 

(2) burning hydrocarbons (benzene and xylene) in air at 
room pressure (samples BE and XY) . The smoke was col- 
lected on quartz substrates. For the UV/vis spectroscopy 
the quartz substrates on which the particles had been col- 
lected were used directly while the dust was scrapped from 
the substrate and embedded in KBr pellets for the IR 
spectroscopy. Bussoletti et al. ( |l987b| ) suggest that the 
extinction efficiencies, Q ox t/a, for the AC samples should 
be corrected by a factor of 5 due to an experimental un- 
derestimation of the pellet density. This cor rectio n gives 
an agreement with t he da ta by Koike et al. ( 1980 ). 

Colangeli et al. (|1995|) measured the extinction effi- 
ciency in the range 40 nm - 2 mm. They produced three 
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different samples; two by arc discharge between amor- 
phous carbon electrodes in Ar and H 2 atmospheres at 10 
mbar (sample ACAR and ACH2 respectively) and one by 
burning benzene in air (sample BE). The samples were de- 
posited onto different substrates for the UV/vis measure- 
ments, while in the IR the samples where prepared both 
on a substrate and by being embedded in KBr/CsI pellets 
and for the FIR measurements the samples were embed- 
ded in polyethylene pellets. These different but overlap- 
ping methods gave the possibility of evaluating the dif- 
ference as a result of embedding the samples in a matrix 
or by having them on a substrate. Colangeli et al. ( 1995 ) 
found that embedding the samples in a matrix introduces 
a systematic error (the matrix effect) while the spectra 
obtained for grains deposited onto a substrate did not suf- 
fer from any matrix effect detectable within the accuracy 
available in the experiment. Therefore the FIR data were 
corrected for the extinction offset introduced by the ma- 
trix. 

Jager et al. ( 1998| ) produced structurally different car- 
bon materials by pyrolizing cellulose materials at different 
temperatures (400° C, 600° C, 800° C and 1000° C), and 
characterised them in great detail. These materials have 
increasing sp 2 /sp 3 ratios making the amorphous carbon 
pyrolysed at 400° C the most " diamond-like" with the 
lowest sp 2 /sp 3 ratio while the amorphous carbon pyro- 
lysed at 1000° C is more "graphite-like" with the highest 
sp 2 /sp 3 ratio. The pyrolysed carbon samples were embed- 
ded in epoxy resin and reflectance of the samples was mea- 
sured in the range 200 nm to 500 /im, making this the first 
consistent laboratory measurement of amorphous carbon 
over the whole spectral range relevant for radiative trans- 
fer calculations of C-rich AGB stars. From the reflectance 
spectra the complex refractive index, m, was derived by 
the Lorentz oscillator method (see e.g. Bohren & Huffman 
1983, Chap. 9). There is a significant difference between 
the two low temperature (400° C and 600° C) and the two 
high temperature samples (800° C and 1000° C). The lat- 
ter two behave very similar to glassy carbon. 

In contrast to grain measurements, the bulk samples 
by Jager et al. ( 1998 ) give the possibility of investiga- 
ting the difference between the influence of the internal 
structure of amorphous material and the morphology of 
the carbon grains. These two properties can be separated 
out due to the careful investigation of the internal struc- 
tures of the four samples and the range of material proper- 
ties that these four amorphous carbon samples span (from 
" diamond- like" to "graphite-like"). 

3. 2. Calculated optical properties of amorphous carbon 

Several authors have used the data of Bussoletti et al. 
( |l987a| ) and Colangeli et al. (1995) to obtain the optical 
constants of amorphous carbon grains. 

Maron ( 1990) us ed the extinction efficiencies of Bus- 
soletti et al. ( 1987a ) (sample AC2) to derive the optical 



constants (n and k) by estimating the complex permit- 
tivity by a combination of the measured absorption ef- 
ficiencies, dispersion formulae and Kramers-Kronig rela- 
tion. The reason for performing these calculations is that 
the optical constants are needed for modelling emission 
properties of grains containing various allotropic carbons 



or having different sizes. Maron (1990) is of the opinion 
that the differences between the primary extinction effi- 



cie ncies o btained by Bussoletti et al. (1987a) and Koike et 



al. ( 1980 ) are real and caused rather by the use of different 
electrodes (amorphous carbon and graphite, respectively) 
than by an underestimation of t he pell et column density 
as suggested by Bussoletti et al. ( 1987b ). Therefore he did 
not int roduce the correction suggested by Bussoletti et al. 
(tL987bD . 

Rouleau & Martin ([199 lh used the AC2 and BE data 
from Bussoletti et al. ( 1987a| ) to produce synthetic optical 



constants (n and k) which satisfy the Kramers-Kronig re- 
lations and highlight the effects of assuming various shape 
distributions and fractal clusters. One of the complica- 
tions in determining these optical properties of amorphous 
carbon material was that in the infrared the extinction 
measurements were done on a sample of sub-micron-sized 
particles, while in the visible and ultraviolet the optical 
constants were obtained by measurements of reflectivity 
and transmission or by electron energy loss spectroscopy 
on bulk samples. These diverse measurements were used 
to produce synthetic optical constants which satisfied the 
Kramers-Kronig relations. 

Preibisch et al. ( 1993 ) used the BE sample from Bus- 
soletti et al. ( 1987a[ ) between 0.1-300 /Ltm and the data of 
Blanco et al. ( 199l| ) between 40-700 /i m, usi ng the same 
technique as used by Rouleau & Martin ( 1991 ) for deriving 
optical constants taking shape and clustering effects into 
account. Preibisch et al. ( 1993 ) extend the available opti- 
cal co nstant s on the basis of the measurements of Blanco 
et al. ( |l99l| ). With these they determine the opacities of 
core-mantle-particles with varying mantle thickness and 
pollution. 



Zubko et al. ( 1996 ) used the ex tinction efficiencies ob- 
tained by Colangeli et al. ( 1995 ) to derive the optical 
constants (n and k) also by use of the Kramers-Kronig 
approach. These data were used to evaluate the possible 
shapes of the amorphous carbon grains in space and the 
possible clustering of the particles. 

In this study w e hav e used the derived optica l 
constants of Mar on (|l990| ), Rouleau & Ma rtin ( |l99l|) , 
Preibisch et al. ( 1993 ), Zubko et al. ( 1996 ) and Jager 



et al. (1199S|) , see Table |j] for details. The extinction ef- 
ficiency data presented in this paper were calculated in 
the Rayleigh approximation for spheres. 

3.3. The nature of silicon carbide 

Thermodynamic equilibrium calcul ation s performed by 
Friedemann (1969a, b) and Gilman ( 1969 ) suggested that 
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Table 1. Comparison of the different laboratory data and a list of authors who have obtained optical constants from 
these data. 



Reference 


Material 


P 


Wavelength 


Designation 


Comments 




name 


(g/cm ) 


interval (fim) 


in this paper 




Bussoletti et al. (1987a) 


AC 2 


1.85 


0.1 - 300 


Bussoletti AC-2 


Arc discharge 


Bussoletti et al. (1987a) 


BE 


1.81 1 


0.2 - 300 


Bussoletti BE 


burning benzene 


Bussoletti et al. (1987a) 


XY 




0.2 - 300 


Bussoletti XY 


burning Xylene 


Colangeli et al. (1995) 


ACAR 


1.87 


0.04 - 2000 


Colangeli AC 


arc discharge in Ar 


Colangeli et al. (1995) 


ACH2 




0.04 - 950 


Colangeli ACH2 


arc discharge in H2 


Colangeli et al. (1995) 


BE 




0.05 - 2000 


Colangeli BE 


burning benzene 


Jager et al. (1998) 


cel400 


1.435 


0.02 - 500 


Jager 400 


most diamond-like 


Jager et al. (1998) 


cel600 


1.670 


0.02 - 500 


Jager 600 




Jager et al. (1998) 


cel800 


1.843 


0.02 - 500 


Jager 800 




Jager et al. (1998) 


cellOOO 


1.988 


0.02 - 500 


Jager 1000 


most graphite-like 


Maron (1990) 


AC2 






Maron 


n & k from Bussoletti et al. (1987a) 


Rouleau & Martin (1991) 


AC 2 






Rouleau 


n & k from Bussoletti et al. (1987a) 


Preibisch et al. (1993) 


BE 






Preibisch 


n & k from Bussoletti et al. (1987a) 


Zubko et al. (1996) 


ACAR 






Zubko 


n & k from Colangeli et al. (1995) 



x Given in Rouleau & Martin (1991). 



SiC particles can form in the mass outflow of C-rich AGB 
stars. The observa tions p erformed by Hackwell ( 1972 ) and 
Treffers & Cohen ( 1974 ) presented the first empirical evi- 
dence for the presence of SiC particles in stellar atmos- 
pheres. A broad infrared emission feature seen in the spec- 
tra of many carbon stars, peaking between 11.0 and 11.5 
fim is therefore attributed to solid SiC particles and SiC is 
believed to be a significant constituent of the dust around 
carbon stars. 

An ultimate proof for the formation of SiC grains in C- 
rich stellar atmospheres was the detection of isotropically 
anomalous SiC grains in primitive meteorites (Bernato- 
wicz et al. 1987 ). Based on isotopic measurements of the 
major and trace elements in the SiC grains and on models 
of stellar nucleosynthesis, it is established that a majority 
of the presolar SiC grains has their origin in the atmos- 
pheres of late-type C-rich stars (Gallino et al. 199C, 1994; 
Hoppe e t al. [1994 ). For re cent reviews see, e.g ., An ders & 
Zinner ( |1993| ), Ott ( |1993| ) and Hoppe & Ott ( |1997| ). 

Detailed laboratory investigations on the infrared 
spectrum of SiC have been presented by the following au- 
thors: Spitzer et al. (1959a, b) performed thin film mea- 
surements on (3- and a-SiC; Stephens ( 198C| ) measu red on 
crystalline /3-SiC smokes; Friedemann et al. ( [1981 ) mea- 
sured two commercially available a-SiC; Borghesi et al. 
( |l985b[ ) investigated three commercially produced a- and 
one commercially produced /3-SiC; Papoular et al. (1998) 
measured two samples of /3-SiC powders, one produced by 
laser pyrolysis and one which was commercially available; 
Mutschke et al. (1999) studied 16 different SiC powders 



which were partly of commercial origin and partly labora- 
tory products (8 cv-SiC and 8 /3-SiC); Speck et al. (|1999[) 



made thi n film measurements of a- and /3-SiC and Ander- 
sen et al. ( 1999 ) have measured the spectrum of meteoritic 
SiC. 



One of the difficulties in interpreting laboratory data 
lies in disentangling the combination of several effects due 
to size, shape, physical state (amorphous or crystalline), 
purity of the sample and possible matrix effects if a ma- 
trix is used. There is a general agreement that grain size 
and grain shape have a crucial influence on the absorption 
feature of SiC particles. This is particularly demonstrated 
by Papoular et al. ( 1998 ) , Andersen e t al. ( 1999] ) and 
Mutschke et al. ( |l999[ ). Papou lar et al. ( |l998| ), Mutschke 



et al. ( 1999| ) and Speck et al. ( 1999 ) have shown that the 
matrix effect does not shift the resonance feature as a 



whole as it was assume d by Friedemann et al. ( 1981 ) and 
Borghesi et al. ( |l985b| ). While Papoular et al. (199S) and 
Mutschke et al. (1999) find that the profile is not shifted 
but altered, Speck et al. ( |1999| ) state that the profile is 
not affected at all, whether a matrix is used in the ex- 
perimental set up or not. The influence of purity of the 
laboratory samples was mainly studied by Mutschke et al. 
( 1999| ) . Another issue considered is the effect of the crystal 
type. Silicon carbide shows pronounced polytypism which 
means that there exist a number of possible crystal types 
differing in only one spatial direction. All these polytypes 
are variants of the same basic structure and can therefore 
be divided into two basic groups: a-SiC (the hexagonal 
polytypes) and /3-SiC (the cubic polytype). It was found 
by Spitzer et al. ( 1959a,b), Papoular et al. (|l998| ), An- 
dersen et al. ( |1999| ) and Mutschke et al. ( |l999| ) that the 
crystal structure of SiC cannot be determined from IR 
spectra, because there is no systematic dependence of the 
ban d profi le on the crystal type. In contrast, Borghesi et 
al. ( 1985bj ) and Speck et al. ( 1999 ) find the contrary result. 



In this paper we have used the average value for bulk 
SiC refl ectan ce spectra of /3-SiC as presented by Mutschke 

= 6.49, wto = 795.4 cm -1 , uj p 
10 to calculate the optical con- 



et al. ( PL999I) with 
= 1423.3 cm" 1 and 7 
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stants n and k, using the one-oscillator model described 
in Mutschke et al. (1999). The damping constant 7 is an 
"ad hoc" parameter, which in a perfect crystal reflects the 
anharmonicity of the potential curve. A damping constant 
of 7 = 10 characterises crystals which are not structurally 
perfect but still far from amorphousness. 

Since there is no systematic dependence of the band 
profile on the crystal type in the data of Mutschke et al. 
( |l999| ), we could just as well have used the data of one 
of their a-SiC samples and would have obtained a similar 
result. 

The optical constants n and k where used to calcu- 
late the extinction efficiency for small spherical grains in 
the Rayleigh limit. Spheres are not necessarily the best 
approximation for the grain shape of SiC particles in C- 
rich AGB stars compared to e.g. a continuous distribution 
of ellipsoids (CDE) as introduced by Bohren & Huffman 
qi983Q . The general appearance of the feature as well as 
the peak position will depend on the grain shape, however, 
common for all grain shapes of SiC are that the feature 
will always fall between the transverse (TO) and the lon- 
gitudinal (LO) optical phonon mode, so the difference will 
be that a spherical grain shape will give rise to a sharper 
and narrower resonance than other grain shape approxi- 
mations. 



4. Dynamical models 

4-1- Modelling method 

To obtain the structure of the stellar atmosphere and cir- 
cumstellar envelope as a function of time we solve the 
coupled system of radiation hydrodynamics and time- 



dependent dust formation (cf. Hofner et al. 1995, Hofner 



& Dorh 1997 and references therein). The gas dynamics is 
described by the equations of continuity, motion and en- 
ergy, and the radiation field by the grey moment equations 
of the radiative transfer equation (including a variable Ed- 
dington factor) . In co ntrast to the models presented in 
Hofner & Dorfi ( 1997 ) we use a Planck mean gas absorp- 
tion coefficient based o n deta iled molecular data as de- 
scribed in Hofner et al. ( 199§| ). Dust formation is treated 



by the so-call ed mo ment method (Gail & Sedlmayr 1988 
Gauger et al 



1990) 



We consider the formation of amor- 
phous carbon in circumstellar envelopes of C-rich AGB 
stars. 

The dynamical calculations start with an initial model 
which represents the full hydrostatic limit case of the grey 
radiation hydrodynamics equations. It is determined by 
the following parameters: luminosity L+, mass M*, effec- 
tive temperature T* and the elemental abundances. We 
assume all elemental abundances to be solar except the 
one of carbon which is specified by an additional parame- 
ter, i.e. the carbon-to-oxygen ratio ec/eo- The stellar pul- 
sation is simulated by a variable boundary (piston) which 
is located beneath the stellar photosphere and is moving 
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Fig. 2. The Planck mean and the Rossland mean calcu- 
lated from the optical constants of amorphous carbon as 
derived from various experiments. See Table |l| for details. 
The Preibisch data would fall right on top of the Rouleau 
data. 



sinusoidally with a velocity amplitude Au p and a period 
P. Since the radiative flux is kept constant at the inner 
boundary throughout the cycle the luminosity there varies 
according to L- m (t) oc Ri n (t) 2 . 

4-- 2. Dust opacities 

The self-consistent modelling of circumstellar dust shells 
requires the knowledge of the extinction efficiency Q ex t of 
the grains, or rather of the quantity Qoxt/a, which is inde- 
pendent of a in the small particle limit which is applicable 
in this context. For the models of long period variables pre- 
sented in Hofner & Dorfi Ql997j ) and Hofner et al. ( |1998| ) a 
fit formula for the Rosseland mean of Q C xt/a derived from 
the optical constants of Maron (1990) was used. 
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Table 2. Comparison of modelling results for different 
dust opacity data: mass loss rate M (in Mo/yr), mean 
velocity at the outer boundary (u) (in km/s), mean de- 
gree of condensation at the outer boundary (/ c ); model 
parameters: L* = 13000 L@, M* = 1.0 Mq, = 2700K, 
ec/£o = 1.4, P = 650 d, A-iip = 4 km/s; 'R' denotes a 
Rosseland mean dust absorption coefficient, 'P' a Planck 
mean; for details see text. 
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model P13C14U4 of Hofner et al. (1998) 



One important point of this paper is to investigate the 
direct influence of Q e xt /a on the structure and wind prop- 
erties of the dynamical models. Therefore we have com- 
puted Rosseland and Planck mean values of Q C xt/a (see 
Fig.||) based on various optical constants derived from 
laboratory experiments (see Sect.|| for details about the 
samples). 

For the dynamical calculations presented here we have 
selected the following data sets (see Table [l| for a detailed 
specification): Jager^OO and JagerlOOO (representing the 
extreme cases), Rouleau (closest to the Maron data used in 
earlier models but extending to wavelengths below 1 /um) 
and Zubko. The data of Preibisch are almost identical to 
the data of Rouleau. 

Figure || demonstrates that for a given data set the 
difference between Planck and Rosseland means is rela- 
tively small. This is due to the fact that amorphous car- 
bon grains have a continuous opacity with a smooth wave- 
length dependence^. 

4-3. Wind properties 

All models discussed here are calculated with the same set 
of stellar parameters, i.e. L* = 13000 Lq, M* = 1.0 Mq, 
T* = 2700 K, e c /e = 1-4, P = 650 d, Au p = 4 km /s, cor - 
responding to model P13C14U4 in Hofner et al. ( |1998| ). 
The only difference between individual models is the 
adopted mean dust opacity. Most models have been cal- 
culated with Rosseland mean dust opacities to allow us a 
direct comparison with earlier models based on Rosseland 
means derived from the Maron data0. 



Wind properties like the mass loss rate M or the time- 
averaged outflow velocity (u) and degree of condensation 
(f c ) are direct results of the dynamical calculations. The 
Rosseland mean models in Tabic ^ (first group) are listed 
in order of decreasing dust extinction efficiency. Both (u) 
and (f c ) change significantly with the dust data used. 
(f c ) increases with decreasing dust extinction efficiency 
while (u) decreases, reflecting a lower optical depth of the 
circumstellar dust shell (see also Sect.|^). The mass loss 
rates seem to show a weak overall trend but it is doubtful 
whether the differences between "neighbouring" models 
in Table || are significant. Since the mass loss rate varies 
strongly with time the average values given in the table 
are more uncertain than the ones for the velocity and the 
degree of condensation which both do not show large vari- 
ations with time. 

The behaviour of the wind properties can be explained 
in the following way: The stellar parameters of the models 
presented here were chosen in such a way that the models 
fall into a domain where dust formation is efficient and the 
outflow can be easily driven by radiation pressure on dust 
(luminous, cool star, relatively high C/O ratio). In this 
case, the mass loss rate is essentially determined by the 
density in the dust formation zone which mainly depends 
on stellar and pulsation parameters (see e.g. Hofner & 
Dorfi 1997). Therefore it is not surprising that the mass 
loss rates of the different models are quite similar. 

On the other hand, in a self-consistent model, the 
degree of condensation (dust-to-gas ratio) depends both 
on the thermodynamical conditions in the region where 
the dust is formed and on the specific grain opacity. The 
higher the mass absorption coefficient the faster the ma- 
terial is pushed out of the zone where efficient dust for- 
mation and grain growth is possible. Therefore the degree 
of condensation decreases with increasing dust absorption 
coefficients (i.e. higher radiative pressure) as grain growth 
is slowed down by dilution of the gas. Note that even 
in the model with the lowest dust absorption coefficient 
(DJ4R) the condensation of "free" carbon (i.e. all carbon 
not locked in CO) is far from complete ((f c ) < !)■ 

For the two extreme cases (JagerlOOO and JagerJ^OO) 
we have also calculated models with Planck mean dust 
opacities. As shown in Table ^| the wind properties of the 
corresponding Planck and Rosseland models (DJ1P/DJ1R 
and DJ4P/DJ4R) are very similar (if the differences are 
significant at all, see above). The two Planck mean mo- 
dels fit nicely into the dust extinction efficiency sequence 
discussed before for the Rosseland mean models. 

As demonstrated in many earlier papers (e.g. Winters 
et al. |1994|; Hofner & Dorfi |1997[) the dust formation in 



1 In contrast, the two means may differ by orders of magni- 
tude for gas opacities in case of molecular line blanketing. 

2 Note however th at al l models use Planck mean gas opacities 
as discussed in Sect. 4.1. 



dynamical models is not necessarily periodic with the pul- 
sation period P. In general the models are multi- or non- 
periodic in the sense that the dust formation cycle is a 
more or less well defined multiple of P. In the models dis- 
cussed here, a new dust shell is formed about every 5-6 
pulsation periods. 
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Fig. 3. Spectral energy distributions of a minimum phase 
model of DROR where the spectra have been calculated 
with different dust data: we compare Rouleau (consistent 
calculation) with JagerlOOO (upper), with Zubko (middle) 
and Jager400 (lower). For details see text. 



In contrast to the hotter inner regions, the atmospheric 
structure below about 1500 K does not repeat after each 
pulsation cycle but is more or less periodic on the dust 
formation time scales which span several pulsation cycles. 
While it is easy to compare the time-averaged wind pro- 
perties of the models it is much more problematic to find 
comparable "snapshots" in different model sequences for 
the discussion of observable properties presented in the 
next section. Therefore in all cases involving detailed ra- 
diative transfer on top of given model structures we have 
decided to show statistical comparisons including several 
maximum and minimum phase models of each sequence. 

5. Spectral energy distributions and synthetic 
colours 

5.1. Frequency- dependent radiative transfer 

The dynamical calculation yields the structure of the at- 
mosphere and circumstellar envelope (density, tempera- 
ture, degree of condensation, etc.) as a function of time. 
The time-independent radiative transfer equation is solved 
for each frequency separately along parallel rays to ob tain 
spectral energy distributions (Windsteig et al. 1997 and 
references therein). The grey gas opacity (Planck mean) 
is taken directly from the dynamical models. The dust 
opacity is calculated from the optical properties of amor- 
phous carbon and in some cases SiC (see Sect. || for the 
description of the different dust data) . 



calculated using the optical constants of amorphous car- 
bon of Rouleau, Zubko, Jdger^OO and JagerlOOO on top of 
a fixed atmospheric structure. Two different kinds of SEDs 
were calculated; (1) fully consistent ones where the same 
amorphous carbon data were used in the dynamical model 
and in the SED calculations and (2) "inconsistent" ones 
where we used different dust opacity data for the detailed 
radiative transfer on top of the same dynamical model 
structure (fixed spatial distribution of density, tempera- 
ture, degree of condensation). The latter spectra enable 
us to distinguish between the effect of the various dust 
data in the radiative transfer calculation and the effect on 
the model structure. 

Figure H shows the result for the SEDs based on a mini- 
mum phase model of the DROR model sequence. The full 
line always denotes the SED of the consistent model where 
the Rouleau data were used for the underlying dynamical 
model as well as for the calculation of the spectrum. In 
the upper panel the (inconsistent) JagerlOOO spectrum 
(dotted) calculated on top of the same Rouleau model is 
shown in addition to the consistent Rouleau spectrum. 
The middle panel shows the same for Zubko and the lower 
panel for Jdger^OO. The effects of the different dust data 
for a given structure compared to a consistent model using 
the Rouleau data can be summarised as follows: 

— JagerlOOO: the spectrum has a lower flux in the short 
wavelength region (0.5 to 5 pan) and the maximum at 
longer wavelengths. The lower flux level in this region 
is due to the fact that Q cx t/a for the JagerlOOO data 
is higher than for the Rouleau data, therefore we have 
a higher total dust opacity which results in less flux 
coming out. The shift of the maximum is also due to 
the higher dust opacity in the JagerlOOO case. 

— Zubko: the spectrum has a lower flux level in the 
short wavelength region and the maximum is shifted 
to longer wavelengths, but not as far as the JagerlOOO 
spectrum. This is due to the fact that JagerlOOO is a 
more "graphite- like" amorphous carbon dust than the 
Zubko material. 

— JdgerJ^OO: the spectrum has a comparable flux all over 
the spectrum and the maximum at slightly shorter 
wavelengths. The slightly higher flux level around 
the maximum results from the lower Qext/a of the 
Jdger^OO data which is due to its more "diamond- like" 
nature compared to the Rouleau data. In the wave- 
length region where the maxima of the spectra lie, the 
two data sets are very similar, therefore the maxima 
of the SEDs do not differ much in wavelength. 

Note that for the "inconsistent" SEDs the total flux 
may differ from the value of the consistent models. 



5.2. Spectral energy distributions 

To investigate how the different dust data influence the re- 
sulting spectral energy distributions (SEDs), spectra were 



5.3. Spectral energy distributions including SiC 

The analysis of mid-IR carbon star spectra indicates that 
SiC is the best candidate to reproduce the observations 
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Table 3. Mean dust optical depths at a few selected wave- 
lengths 
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Fig. 4. Spectral energy distributions of a minimum phase 
model of DJ1R where the spectra have been calculated 
with amorphous carbon alone (dotted), with a ratio 
X am c:Xsic of 9:1 (dashed) and a ratio of X amC :X S ic of 
4:1 (solid). 



around the 11 /im region. We have therefore considered 
SiC as an additional dust component (see Sect. 3.2 for de- 
tails). The formation of SiC is not included in the self- 
consistent model calculations because (1) little is known 
about the condensation process and (2) because we do 
not expect that SiC will have a significant influence on 
the model structures. We use either Planck or Rosse- 
land means for the model computations and SiC will con- 
tribute only in a very narrow wavelength region with small 
amounts to these mean opacities compared to amorphous 
carbon. 

The effect of SiC as dust component is described in a 
qualitative manner. The dust opacity k& for each wave- 
length is calculated from 



;-y a mC+KSiC^SiC 



where X; are the fractional parts of amorphous carbon 
and SiC, respectively, and where K a mC an d ^sic are the 
opacities of carbon and SiC. Figure || shows how a mixture 
of dust grains consisting of amorphous carbon (JdgerlOOO 
data) and SiC modifies the SED around 11.3 /im. Two 
different ratios of X am c : Xs;c, 4:1 and 9:1 were adopted. 
The higher the amount of SiC, the stronger is the 11.3 /xm 
feature (see inset of Fig. |J) . Another choice of grain shape 
than spherical for the SiC particles, would result in a 
broader and weaker feature. 



5.4- Synthetic colours 

For a comparison of the consistent spectra (model struc- 
ture and spectra computed with the same dust data) we 
have calculated synthetic J, H and L colours as well as 
the IRAS 12 fim colour. In a (J— H) versus (H— L) dia- 
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gram (Fig. ||a) the models based on different amorphous 
dust data fall into distinct regions. 

The models calculated with the JdgerlOOO dust data 
have the reddest colours in (H— L). The Jdger 400 colours 
are the bluest, while Rouleau and Zubko lie in between. 
In (J— H) models with the Zubko data are reddest and 
the others do not differ much. The reason for the diffe- 
rent slopes of Rouleau and Zubko compared to both of the 
Jdger data sets is that in these cases the maxima of the 
SEDs are changing between the J and the H filter depen- 
ding on the phase. The maxima of the SEDs resulting from 
the JdgerlOOO model structures are always at longer wave- 
lengths and the ones of the Jdger 4-00 structures lie mainly 
in one filter. From Fig. |^b, which shows the "inconsistent" 
colours based on model DROR (structure was calculated 
with the Rouleau data and the spectra with other dust 
data) in addition to the consistent Rouleau colours, it is 
clear that the influence of the different dust data used in 
the radiative transfer calculation is much stronger than 
the effect of the underlying hydrodynamic model struc- 
ture. 

Note that in Fig.^ only maximum and minimum 
phases are shown. Other phases would fill in the gaps be- 
tween successive extremes. The colours are strongly re- 
lated to the formation of a new dust shell which takes 
place every 5 to 6 pulsation cycles (see Sect. 4.3). After 



this time scale the colours match very closely the ones 
of the preceding dust formation cycle as shown in Fig. ^| 
When connecting the succeeding points it can be seen that 
they form a spiral. The minima (circles) are always redder 
than the following maxima (asterisks). 

To investigate also the mid-IR properties of the models 
we calculated the 12 /im colour. A (L— [12]) vs. (H— L) dia- 
gram (Fig.||c) shows, that again the consistent colours fall 
into distinct regions. The sequence in (L— [12]) {Zubko - 
Rouleau - JdgerlOOO - Jdger 4-00) is a sequence of decreas- 
ing optical depths. Table [| lists the mean dust optical 
depths for a few selected wavelengths. In Fig.|^d the in- 
consistent colours based on the model structure of DROR 
are shown for comparison. 

From Fig.[| we can infer that the influence of the dif- 
ferent amorphous carbon dust data used in the radiative 
transfer calculation is much stronger than the effect of 
the model structures. The colours resulting from the same 
dust data fall approximately into the same region of a two- 
colour-diagram, whether they are calculated on top of the 
corresponding model structure (upper panel) or a fixed 
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This plot shows that the influence of the different dust data used in the radiative transfer calculation is much stronger 
than the effect of the model structures. 




Fig. 6. Upper panel: (J— H) vs. (H— L) for maxima (as- 
terisks) and minima (circles) of the DJ1R model. When 
connecting the points following each other one can see that 
the result is a spiral, the minima are always redder than 
the following maxima. The colours are strongly related to 
the formation of a new dust shell which takes place every 
5 to 6 pulsation cycles. After this time scale the colours 
match very closely the ones during the formation of the 
prior dust shell. Lower panel: Same as above, only (H— L) 
vs. (L— [12]) is shown. 



model sequence (lower panel). This applies for all data 
sets. 

In addition we computed synthetic colours for the 
Maron and the Preibisch data, but they are not shown in 
Fig. U to avoid overlaps in the plot. The Preibisch colours 
fall into the same region as the Zubko colours as one would 
expect because the Q cyi t/a are very similar. The same ap- 
plies for the Maron and the Rouleau data. Only the J flux 
differs between Maron and Rouleau, the reason being that 
the Maron data do not extend below 1 /an and therefore 
the corresponding contribution to the J filter is lacking. 

6. Summary and conclusions 

Carbon bearing grains are expected to form in the out- 
flows of C-rich AGB stars. The two most common types 
of carbon grains in these stars are expected to be amor- 
phous carbon and SiC grains. We have investigated the 
direct influence of different dust optical properties on the 
wind characteristics and the resulting observable proper- 
ties of the dynamical models. 

The term amorphous carbon covers a wide variety 
of material properties from "diamond-like " to "graphite- 
like" . We ha ve us ed n and k data of M aron (1990 ), Rouleau 
& Martin ( |l99l| ), Preibisch et al. fll993| ), Zubko et al. 
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( |1996| ) and Jager et al. (1998) to investigate the influence 
of different types of amorphous carbon, on the structure 
and the wind properties of dynamical models. The Rosse- 
land and Planck mean values of Q ex t / & used in the model 
computations were calculated from the Rayleigh approxi- 
mation for spheres. The difference between the Planck and 
the Rosscland means is relatively small for the amorphous 
carbon data because the grains have a continuous opacity 
with a smooth wavelength dependence. 

In our test models, both the outflow velocity (u) and 
degree of condensation (/ c ) change significantly with the 
dust data used. (/ c ) increases with decreasing dust ex- 
tinction efficiency while (u) decreases, reflecting a lower 
optical depth of the circumstellar dust shell. The mass 
loss rate is, however, not significantly influenced by the 
use of different dust data. 

On top of the structures resulting from the dynamic 
calculations we have performed detailed radiative transfer 
calculations to obtain the spectral energy distribution of 
the circumstellar dust shells. Regarding infrared colours, 
the influence of the different dust data used in the radia- 
tive transfer calculation is much stronger than the effect of 
the underlying hydrodynamic model structure. However, 
this should not be used as an excuse for fitting observa- 
tions by arbitrarily choosing the optical properties of the 
dust grains for a given model structure. In a consistent 
model the dynamical properties (e.g. outflow velocities) 
and the optical appearance of the circumstellar envelope 
are related in a complex way. 

The influence of including SiC grains is that the 
11.3 /xm feature appears in the spectral energy distribu- 
tion of the models. How much SiC should be "mixed" into 
a model to repro duce the 11.3 /im feature observed (e.g. 
class 4 in Goebel 1995 ) will very much depend on the as- 
sumptions which are made about the size and shape of 
the S iC grains which enter into the model ( Papou lar et al. 

Mutschke et al. 



199? ; Andersen et al. 1999 



1999) 
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